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Abstract 

One of the most important non-dimensional parameters 
used in ice-accretion modeling and scaling studies is the 
freezing fraction defined by the heat-balance analysis of 
Messinger. For fifty years this parameter has been used 
to indicate how rapidly freezing takes place when su- 
per-cooled water strikes a solid body. The value ranges 
from 0 (no freezing) to 1 (water freezes immediately on 
impact), and the magnitude has been shown to play a 
major role in determining the physical appearance of 
the accreted ice. Because of its importance to ice 
shape, this parameter and the physics underlying the 
expressions used to calculate it have been questioned 
from time to time. Until now, there has been no strong 
evidence either validating or casting doubt on the cur- 
rent expressions. 

This paper presents experimental measurements of the 
leading-edge thickness of a number of ice shapes for a 
variety of test conditions with nominal freezing frac- 
tions from 0.3 to 1.0. From these thickness measure- 
ments, experimental freezing fractions were calculated 
and compared with values found from the Messinger 
analysis as applied by Ruff. Within the experimental 
uncertainty of measuring the leading-edge thickness, 
agreement of the experimental and analytical freezing 
fraction was very good. It is also shown that values of 
analytical freezing fraction were entirely consistent 
with observed ice shapes at and near rime conditions: 
At an analytical freezing fraction of unity, experimental 
ice shapes displayed the classic rime shape, while for 
conditions producing analytical freezing fractions 
slightly lower than unity, glaze feahires started to ap- 
pear. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 


c p Specific heat of air, cal/g K 
c P ,ws Specific heat of water at the surface tempera- 
ture, cal/g K 

CL Tunnel centerline (mid-span) position 
d Cylinder diameter or twice the leading-edge 

radius of airfoil, cm 

D v Diffusivity of water vapor, cm 2 /sec 

h c Convective heat-transfer coefficient, 

cal/s m 2 K 

h G Gas-phase mass-transfer coefficient, g/sm 2 

k a Thermal conductivity of air, cal/s mK 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

LWC Cloud liquid-water content, g/m 3 

MVD Water droplet median volume diameter, pm 
n Freezing fraction, dimensionless 

n a Freezing fraction calculated using Messinger 

analysis, dimensionless 

n e Freezing fraction from leading-edge ice thick- 
ness, dimensionless 

Nu Nusselt number, dimensionless 

p Static pressure, Nt/m 2 

p w Vapor pressure of water in atmosphere, Nt/m 2 

p m Vapor pressure of water at the icing surface, 

Nt/m 2 

r Recovery factor, dimensionless 

Re Reynolds number of model, dimensionless 

Res Reynolds number of water droplet, dimen- 

sionless 

Sc Schmidt number, dimensionless 

SLD Super-cooled large droplet 

tf Freezing temperature of water, °C 

t s Surface temperature, °C 

t Temperature, °C 

T Absolute temperature, K 

V Free-stream velocity of air, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

We c Weber number based on model size and air 
properties, dimensionless 

Wei, Weber number based on water-film thickness 
and water properties, dimensionless 
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pa Collection efficiency at stagnation line, dimen- 
sionless 

A Ice thickness at stagnation line, cm 

{ i> Water droplet energy transfer parameter, °C 

A Water droplet range, m 

A Sl okes Water droplet range if Stokes Law applies, m 

Af Latent heat of freezing of water, cal/g 

A v Latent heat of evaporation of water, cal/g 

fi Viscosity of air, g/m s 

0 Air energy transfer parameter, °C 

p Air density, g/m 3 

Pi Ice density, g/m 3 

p w Liquid water density, g/m 3 

r Accretion time, min 

Subscripts 

st Static 

tot Total 


Introduction 


One of the most important non-dimensional parameters 
used in ice-accretion modeling and scaling studies is the 
freezing fraction, n. Since Messinger 1 introduced the 
term in his surface heat balance analysis fifty years ago, 
this parameter has been shown to have a strong effect 
on the shape of ice accretions. 2, 3 It is important that 
analytical expressions used to calculate it include all 
significant physical processes. 

Messinger defined the freezing fraction as the fraction 
of water flux entering a control volume that freezes 
within the control volume. Along the stagnation line of 
a cylinder or airfoil water reaches the surface only by 
impingement. The portion that freezes is represented 
by the accumulated ice thickness. Thus, the non- 
dimensional leading-edge ice thickness, A, is: 

4 = nA *A (!) 

a 


Here, d is the cylinder diameter or twice the leading- 
edge radius for an airfoil. For the NACA 0012, the 
leading-edge radius is 0.0158c, where c is the chord. A c 
is the accumulation parameter 


LWC Vx 
d p t 


(2) 


and /3 0 is the collection efficiency of water droplets at 
the stagnation line. The Langmuir and Blodgett 4 
method for calculating /3 0 will be presented below. In 
this paper we shall use n e to refer to freezing fractions 
calculated from experimental ice shapes using equation 
(1) and n a for analytical freezing fractions calculated 
from the equations given below. 


For scaling studies, only the stagnation value of the 
freezing fraction is matched between scale and refer- 
ence conditions. It is assumed that if a match is 
achieved there, the freezing fraction will match every- 
where over the bodies of interest. The value of this 
stagnation freezing fraction as determined from analysis 
must be able to pass two tests, (i) It must lead to values 
consistent with the experimental ice thickness at the 
stagnation line (from eq. (1)), and (ii) it must also lead 
to values consistent with the ice shape at and near the 
leading edge. Test (ii) means that for an analytical 
freezing fraction of 1, ice shapes must be rime, and if n a 
is calculated to be less than 1, ice shapes must show 
glaze features at the stagnation line. If either of the 
requirements (i) or (ii) is not met, the analysis must be 
deficient. 

In 1988 Bilan in' analyzed Ruff 2 ice tracings and calcu- 
lated experimental freezing fractions using equation (1). 
Only for high freezing fractions did he find reasonable 
agreement between n e and the Messinger-analysis 
freezing fraction, n a , and he concluded that at low n the 
currently-used expression may not adequately account 
for water-film effects. However, the value of n e found 
by equation (1) is strongly dependent on the uncertain- 
ties inherent in recording and measuring small ice 
thicknesses and also somewhat affected by how accu- 
rately the L WC is known for the icing tests. These po- 
tential problems were not addressed by Bilanin who 
worked with a limited data set. 


This paper presents experimental measurements of the 
leading-edge thickness of a number of ice shapes re- 
corded in the NASA Glenn Icing Research Tunnel 
(IRT) over a range of freezing fractions from 0.3 to 1.0. 
From the thickness measurements, n e was found from 
equation (1) and compared with the value of n a found 
using the method of Ruff, which was based on the 
Messinger surface heat-balance analysis. 

Analysis 

Stagnation Collection Efficiency 

Each of the parameters of importance to this study, in- 
cluding freezing fraction, will be used in the form pre- 
sented by Ruff 2 and subsequently applied by Anderson ’ 
in recent scaling studies. 


The stagnation collection efficiency, j3 0 , in equation (1) 
was calculated from the expression given by Langmuir 
and Blodgett: 


1.40(AT o - i/s) 84 
1 + L40(7V 0 -1/8) 84 


( 3 ) 


where the modified inertia parameter, K 0 , is: 
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(4) 


Here A/A Sto k es is the droplet range parameter, defined as 
the ratio of actual droplet range to that if Stokes drag 
law for solid-spheres applied. It is a function only of 
the droplet Reynolds number. Res and was tabulated by 
Langmuir and Blodgett. The following fit to the Lang- 
muir and Blodgett tabulation was used in this study: 


A 


f 0.8388 + 0.0014837^" 


+ 0 . 1847^7 


where the droplet Reynolds number is 
V MVD p w 


Re s =- 


P 


The inertia parameter, K, in equation (4) is 
Pw MVD 1 V 


K = 


18 d f.i 


(5) 


( 6 ) 


(7) 


Analytical Stagnation Freezing Fraction 

The freezing fraction, from Messinger’s surface energy 
balance, can be written in the form, 


n a 


p.ws 





( 8 ) 


difference in the values resulting from the two. Equa- 
tion (10) was used in this study. 

The second term in equation (10) accounts for evapora- 
tion. The gas-phase mass-transfer coefficient, /? G , can 
be found from 






where the Prandtl number is 



and the Schmidt number is given by 


Sc 


P 

pD v 


( 12 ) 


(13) 


(14) 


The convective heat-transfer coefficient, h c , was deter- 
mined from the Nusselt number: 


h c = —Nu (15) 

d 

The evaluation for Nu at the stagnation line of a cylin- 
der with laminar flow was reported by Kreith: x 

Nu=l.l4Pr 0A Re 0 - 5 (16) 

with 


This expression does not include heat lost from the sur- 
face due to water runback or to conduction into the 
model surface. The latter component of the heat bal- 
ance is not likely to be important except when accretion 
first starts or soon after. The individual terms in this 
expression are (/). the water energy transfer parameter; 
ft the air energy transfer parameter; and b, the relative 
heat factor, introduced by Tribus, et. al. 6 These pa- 
rameters are defined as 


<P = t f -t sl 


2 c n 


(9) 


e = \ 


t. -L, 


rV 1 


2c 


P J 


Pww _ Ptot Pw 
T S t T lot Pst 
1 Plot _ Pwv[ 
0.622 T tot T st j 


A (10) 


Re = 


Vdp 


P 


(17) 


The air properties used in equation (16) should be based 
on the film temperature. For temperatures of interest to 
icing, Pr is effectively constant, and equation (16) can 
be simplified to 

Nu = 0.992 fle 05 (18) 

Poinsatte 9 measured heat-transfer coefficients at the 
stagnation line of an NACA 0012 airfoil in the IRT. 
Poinsatte defined his parameters using the chord as the 
length, but to be consistent with the practice in this 
study, they can be rewritten in terms of d, using d = 
0.0316c for the NACA 0012 airfoil. Poinsatte’s expres- 
sion for Nu at 0° angle of attack is then 

Nu = l.lORe 0472 (19) 


b = 


LWC V /? 0 c p 

K 


(ii) 


and at 4° AOA, 

Nu = 0.903 Re 0498 


( 20 ) 


Equation (10) gives the form of ft employed by Ruff, 
which includes compressibility effects. Chaipin and 
Fasso 7 and others have used a simpler form that ignores 
compressibility, but for icing conditions there is little 


These correlations are remarkably close to the expres- 
sion for cylinders in equation (18), supporting the ap- 
proach used here of applying cylinder expressions to 
airfoils with the airfoil leading-edge radius substituted 
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for cylinder radius. Equation (16) was used in the pre- 
sent study. 

IRT and Test Description 

The icing test results used in this investigation were 
performed in the NASA Glenn IRT. The IRT is a 
closed-loop, refrigerated, sea level tunnel with a test 
section of 1.8 by 2.7 m. The IRT has 10 spray bars. 
The cloud from these spray bars is calibrated periodi- 
cally; for the tests reported here the calibration com- 
pleted in the summer of 2000 10 was used. 

The models used for these tests were NACA 0012 61- 
cm-span aluminum airfoil sections with chords of 80.0, 
53.3, 35.6 and 26.7 cm. They were mounted vertically 
between splitter plates at the center of the IRT test sec- 
tion as shown typically in figure 1. Guide lines for lo- 
cating ice-tracing templates were marked at the leading 
edge of each model at tunnel center, ±2.5 cm and ±5 cm 
from the center. All tests were run at 0° AOA. 

In preparing for a test, the temperature and airspeed in 
the test section and the air and water pressures on the 
spray manifolds were set. When these conditions stabi- 
lized at the set points, the spray nozzle valves were 
opened to initiate the spray. The IRT spray system 
reaches steady state in a few seconds, eliminating the 
need to shield the models during spray stabilization. At 
the completion of the spray duration, the nozzle valves 
closed and the tunnel fan was brought to a full stop. 
Personnel then entered the test section to record the ice 
shape through tracings and photographs. 

For tracings, a heated ice knife cut from 0.16- to 0.32- 
cm-thick stainless-steel sheet was inserted into the lead- 
ing edge of the ice to melt a thin horizontal slit down to 
the model surface. A cardboard template cut out to 
conform with the shape of the model near the leading 
edge was placed into this gap and a cross section of the 
ice shape traced onto the template with a pencil. Trac- 
ings were taken at two locations: the vertical center of 
the tunnel (91 cm from the floor) and at 2.5 cm above 
the center. The tracings were digitized and the x-y co- 
ordinates for each ice shape recorded. The coordinates 
were normalized by the model chord when preparing 
presentation plots. 

Ice shapes used in this study were obtained in scaling 
studies made in March 2001 and March and October 
2002. Each of those test series included SLD and Ap- 
pendix-C test conditions, and the 2001 scaling results 
have been reported previously." 12 To limit the analy- 
sis, only the Appendix-C subset of those tests was used 
for the present work. Three rime tests were eliminated 
from that subset because the recorded ice shapes had 
been obtained at conditions for which erosion had been 
seen in the past. What remained were forty sets of ice 


shapes at conditions producing a range of analytical 
freezing fractions from 0.28 to 1. 

The leading-edge ice thickness for each ice shape was 
determined by examining the x-y coordinate data near 
the y = 0 location. Because the tests were performed at 
0° AOA and the origin of the coordinate system was at 
the leading edge of the clean model, the absolute value 
ofx at y = 0 would ideally represent the ice thickness at 
the stagnation line. However, slight misalignments of 
the model with the flow sometimes displaced the actual 
stagnation line slightly above or below y = 0. 
Therefore, the local minimum (for glaze ice) or 
maximum (for rime) absolute value ofx near y = 0 was 
selected as the ice thickness. Ice shape differences be- 
tween the two tracing locations typically did not appear 
to be significant, but differences in ice thickness at the 
leading edge were usually measurable. 

Uncertainty Analysis 

Tunnel and cloud conditions were recorded every 1, 3 
or 5 sec (depending on spray duration) over the length 
of the spray. Reported conditions are the time averages 
of these data histories. Estimates of the uncertainty in 
the reported average conditions were made by consider- 
ing fluctuations of the values over time, possible in- 
strument errors, uncertainties in tunnel calibration of 
MVD and L WC and differences in measurements from 
one location to another in the test section. 

Total temperature was measured by 24 thermocouples 



Figure 1 . 61-cm-Span Model Installed in IRT Test 
Section. 
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distributed over the turning vanes just upstream of the 
spray bars. These values were then averaged at each 
time interval to determine a temperature-time history. 
Including inherent uncertainty in the thermocouples 
themselves, it was estimated that the total uncertainty in 
temperature was on the order of ±0.5°C. 

Test-section velocity was determined from the average 
total and average static pressures from 2 pitot-static 
probes, one located on the north wall and one on the 
south wall of the tunnel at the entrance to the test sec- 
tion. The maximum computed difference in velocities 
indicated by the two probes was ±0.25%. Pressures 
also fluctuated modestly with time such that the varia- 
tion in velocity over the spray duration was less than 
3%. 

The most significant contribution to L WC uncertainty is 
the calibration process, and scatter in the calibration 
data suggests this uncertainty was about ±12%. 

Most of the uncertainty in drop-size measurement came 
from interpreting results from the instruments used in 
calibration. The net uncertainty in MVD was estimated 
at ±12%. This uncertainty applies only to the particular 
instruments used (OAP and FSSP probes) in the 1RT 
calibration and do not suggest a reference to an absolute 
value, which is unknown. 

These estimated uncertainties in the test conditions pro- 
duced approximate uncertainties of ±2% in j3 0 , ±12% in 
A c and ±15% in n. 

Results and Discussion 
Leading-Edge Thickness 

The measured leading-edge thickness values, corre- 
sponding freezing fractions using equation (1) and the 
conditions for each test are listed in Table I in chrono- 
logical order. Thickness values and n e for both mid- 
span (CL) and 2.5 cm above mid-span (CL ± 2.5 cm) 
are given. The absolute difference between the 
CL±2.5-cm and the CL thickness for each test varied 
from 0 to 0.14 cm. Also included in the table are the 
values of / 3 0 , A c and n a computed using equations (2) to 
(17). 

Sample ice tracings are reproduced in figure 2. The 
mid-span (CL) shape is shown with a shaded cross- 
section while the 2.5-cm-above mid-span (CL±2.5 cm) 
is represented with a solid line. The spray conditions 
for figure 2 (a) produced a glaze accretion with an n a of 
0.28. While the tracings for the two locations looked 
very much alike, the leading-edge thicknesses for the 
two differed by 18.8%. For low freezing fractions the 
features of the main ice shape tended to be somewhat 
random in nature; therefore, the leading-edge location 
and thickness were often difficult to define. Further- 


more, because the ice thickness was small for low 
freezing fractions, any inherent inaccuracies in the trac- 
ing process would have tended to produce larger per- 
centage differences from tracing to tracing than would 
occur at higher freezing fractions. For the tests of Ta- 
ble I, when n a was less than 0.3 the CL±2.5-cm thick- 
ness varied from -15.8 to 30.2% of the CL value. The 
highest relative difference occurred for a test using the 
smallest model, for which the inherent inaccuracies of 
manually tracing the ice were the greatest. 

Figure 2 (b) gives the ice shapes recorded for a glaze 
accretion with a freezing fraction of 0.52. The leading- 
edge features were better defined than for the shapes 
shown in figure 2 (a), and it was relatively easy to pick 
out the leading edge thickness from the ice-shape coor- 
dinates. Here the ice tracings for the two locations 
yielded thickness values that differed by 10.2%. 



x/c 

(a) Ice Shapes from 3-23-01 Run 2. Glaze Ice with n a 


= 0.28. 



x/c 

(b) Ice Shapes from 3-23-01 Run 7. Glaze Ice with n a 
= 0.52. 

Figure 2. Sample Ice Shapes. See Table I for Test 
Conditions. 
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Table I. 


Test Conditions, Leading Edge Thicknesses and Freezing Fractions 


Date/Run 

C, 

cm 

tsh 

°c 

V, 

m/s 

MVD , 

pm 

LWC, 

g/m 3 

% 

min 

Po 

A c 

n* 

4 

CL 

cm 

CL + 
2.5 cm 

n 

CL 

e 

CL + 
2.5 cm 

3-13-01/5 

80.0 

-7.1 

57.5 

45.9 

0.915 

13.8 

0.827 

p.879 

0.292 

1.21 

1.26 

0.308 

0.322 

3-13-01/7 

80.0 

-7.0 

48.5 

43.2 

0.979 

Pi 4.5 

O.SOT 

P.780 

0.310 

1.17 

pL14 

0.325 

0.315 

3-14-01/5 

80.0 

-9.9 

57.7 

44.2 

0.899 

P 13.9 

0.820 

p.868 

0.42 1~] 

1.63 

1.58 

0.420 

0.409 

3-14-01/6 

80.0 

-9.8 

47.6 

41.4 

0.958 

Pi 4.9 

0.789 

p.757 

0.453 

1.52 

1.58 

0.433 

0.45 P] 

3-14-01/10 

80.0 

-12.6 

58.1 

43.4 

0.889 

14.0 

0.817 

1.871 

0.541 

1.84 

1.98 

0.477 

0.513 

3-15-01/4 

80.0 

P-14.9 

57.8 

43.2 

0.891 

Pi 4.0 

0.816 

M.865 

0.637 

2.32 

2.26 

0.604 

0.588 

3-15-01/10 

80.0 

-18.7 

57.9 

43.6 

0.885 

Pi 4.1 

0.817 

P.872 

0.791J 

2.87 

3.04 

0.74lP 

0.787 

3-19-01/5 

53.3 

-7.2 

66.8 

39.9 

0.996 

7.3 

0.857 

P.885 

0.284 

0.70 

0.59 

0.255 

0.215 

3-19-01/6 

53.3 

-9.9 

66.8 

39.9 

0.996 

7.3 

0.857 

P.884 

0.403 

0.98 

0.95 

0.359 

0.35H 

3-20-01/1 

53.3 

-7.1 

66.8 

40.0 

0.997 

7.3 

0.858 

P.886 

0.279 

0.6T1 

0.69 

0.224 

0.253 

3-20-01/5 

53.3 

-12.7 

66.9 

40.0 

0.995 

7.3 

0.858 

P.887 

0.517 

pL13 

1.2 lj 

0.416 

0.444 

3-20-01/8 

53.3 

-15.0 

66.9 

40.0 

0.995 

7.3 

0.858 

P.886 

0.611J 

1.37 

1.33 

0.503 

0.486 

3-20-01/9 

53.3 

-18.9 

67.0 

40.0 

0.995 

7.3 

0.858 

P.887 

0.766 

2.16 

2.17 

0.792 

0.797 

3-21-01/4 

53.3 

r- 26 . n 

66.8 

40.0 

0.997 

7.3 

0.858 

P.886 

P .000 

2,61 

2.68 

0.957 

0.984 

3-22-01/4 

26.7 

P-19.6 

86.fi 

23.6 

P”fl82 

2.4 

0.859 

P.896 

0.775 

0.96 

1.03 

0.699 

0.752 

3-22-01/5 

26.7 

P-15.4 

85.8 

23.6 

1.173 

2.4 

0.858 

p.876 

0.610 

0.8 G 

0.80 

0.597 

0.589 

3-22-01/7 

26.7 

-13.1 

86.1 

23.5 

1.159 

2.4 

0.858 

1.860 

0.517 

0.67 

0.68 

0.496 

0.507 

3-22-01/8 

26.7 

P-10.5 

87.4 

23.5 

L124 

2.5 

0.859 

p.907 

0.408 

0.54 

0.56 

0.393 

0.408 

3-22-01/Ifl 

26.7 

-7.7 

88.4 

23.3 

1.065 

2.6 

0.858 

P.900 

0.285 

0.45 

0.42 

0.330 

0.307 

3-23-01/1 

35.6 

-7.7 

87.8 

28.1 

0.930 

4.0 

0.859 

1.901 

0.287 

0.52 

0.63 

0.285 

0.345 

3-23-01/2 

35.6 

-7.3 

77.6 

29.2 

1.062 

3.9 

0.858 

1.871 

0.278 

0.51 

0.60 

0.282 

0.335 

3-23-01/3 

35.6 

P-10.3 

90.2 

27.8 

0.941 

3.8 

0.859 

P.880 

0.395 

o.tTI 

0.72 

0.39f1 

0.398 

3-23-01/4 

35.6 

P-10.2 

78.0 

29.5 

1.080 

3.9 

0.859 

P.912 

0.403 

0.74 

0.78 

0.403 

0.425 

3-23-01/5 

35.6 

P-13.2 

91.4 

27.6 

0.950 

3.7 

0.859 

1.871 

0.516 

0.88 

0.88 

0.487 

0.489 

3-23-01/6 

35.6 

P-12.9 

77.9 

29.3 

1.091 

3.8 

0.859 

P.880 

0.519 

0.87 

0.90 

0.479 

0.499 

3-23-01/7 

35.6 

P-13.0 

77.8 

29.3 

1.092 

3.8 

0.859 

P.879 

0.522 

0.9T1 

1.00 

0.503 

0.554 

3-23-01/8 

35.6 

P-15.5 

92.6 

27.5 

0.952 

3.7 

0.859 

P.899 

0.610 

1.07 

1.06 

0.581 

0.576 

3-23-01/9 

35.6 

-15.1 

77.7 

29.5 

1.098 

3.8 

0.859 

p.888 

0.607 

Pup 

pf 12 

0.603 

0.612 

3-23-01/10 

35.6 

-19.3 

93.3 

27.3 

0.959 

3.6 

0.859 

1.876 

0.760 

1.38 

1.24 

0.764 

0.687 

3-23-01/11 

35.6 

-19.0 

77.2 

29.4 

L108 

3.8 

0.859 

p.894 

0.759 

1.30 

1.36 

0.709 

0.742 

3-23-01/13 

35.6 

P-26.2 

77. 1_\ 

29.4 

1,115 

3.8 

0.859 

P.903 

P.OOO 

1.67 

1.65 

0.910 

0.90 M 

3-12-02/1 

26.7 

-10.2 

P 14.5 

21.4 

1.238 

1.9 

0.862 

2.091J 

0.275 

0.29 

0.38 

0.190 

0.248 

3-12-02/2 

26.7 

-13.2 

[1 14.4 

21.5 

L193 

2.0 

0.863 

2.119 

0.395 

0.49 

0.50 

0.319 

0.325 

3-12-02/6 

26.7 

-16.0 

p 14.6 

21.5 

L169 

2.1 

0.863 

2.184 

0.504 

0.68 

0.61J 

0.426 

0.383 

3-12-02/8 

26.7 

-18.4 

p 14.8 

21.5 

L155 

2.1 

0.863 

2.161] 

0.598 

0.88 

0.82 

0.558 

0.52 M 

10-11-02/8 

53.3 

P-is.n 

51.5 

40.1 

1.129 

PlO.8 

0.84 ip 

2.439 

0.777 

2.75 

2.39 

0.795 

0.69H 

10-11-02/9 

53.3 

P-19.4 

51.5 

40.1 

L129 

PlO.8 

0.84 Ip 

2.439 

0.824 

2.69 

2.87 

0.777 

0.829 

10-12-02/7 

53.3 

P-20.4 

51.4 

40.1 

L131 

PlO.8 

0.840 

2.437 

0.864 

2.77 

2,91 

0.802 

0.844 

10-12-02/8 

53.3 

P-21.7 

51.4 

40.1 

L131 

PlO.8 

0.84 )P 

2.437 

0.914 

3.05 

3.17 

0.883 

0.918 

10-12-02/11 

53.3 

P-25.0 

50.7 

40.2 

1141 

PlO.8 

0.84 P| 

2.429 

p.ooo 

3.52 

3.68 

Pl.02f1 

Pf070 


* n a Calculated Using Nu from equation (16) 


Rime shapes ( n a = 1.0) reproduce the collection- 
efficiency curve around the airfoil because all water 
impinging at a given location will freeze there. Typical 
rime shapes used in this study are given in figures 2 (c). 
The relatively smooth leading edge with a well-defined 
maximum makes it easy to establish the ice thickness 


from a table of the ice coordinates. The two tracings 
compared in the figure had ice thicknesses that differed 
by only 2.8%. 

Rime tests can be used to help validate tunnel LWC 
calibrations. The temperature used for the test of figure 
2 (c) was calculated as the maximum for which fully 
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0.08 



xlc 

(c) Ice Shapes from 3-21-01 Run 4. Rime Ice with n a = 

1 . 00 . 

Figure 2. concluded. 

rime ice would form for the specified model size, and 
test conditions. If the true LWC had been significantly 
higher than indicated from the calibration, evidence of 
glaze ice at the leading edge would have been seen. On 
the other hand, if the true L WC was lower than the cali- 
brated value fully rime accretions would be seen for 
higher temperatures. Thus, for a given set of test condi- 
tions, if there is confidence in the temperature, velocity 
and droplet MVD, a series of tests in which the 
temperature is incrementally increased in small steps 
above the calculated maximum value for rime should 
either validate or identify problems with the L WC cali- 

iKfr 'temperatures higher than the maximum temperature 
for rime, n a will be less than unity. Therefore, glaze 
features should occur at the leading edge for there to be 
consistency between the analytical freezing fraction and 
experimental ice shapes. Results of tests with tempera- 
tures incremented above the maximum for rime will be 
presented later. 

Comparison of Analytical and Experimental Freezing 
Fraction 

The experimentally-derived freezing fractions are plot- 
ted against the Messinger analytical values in figure 3. 
Open symbols were measured from the ice tracings at 
mid-span (CL) and the solid from the tracings 2.5 cm 
above mid-span (CL+2.5 cm). The solid line represents 
perfect agreement of the two freezing fractions and the 
shaded band indicates the limits for ±10% agreement. 
A linear fit to the data is also shown as a dashed line. 

As noted above, there can be significant uncertainty in 
the ice thickness values found from tracings at low 
freezing fractions. The data scatter in figure 3 at the 
lowest freezing fraction tested was a consequence of 
this uncertainty. The linear fit to the data fell within the 



Freezing Fraction from Analysis. n a 

Figure 3. Experimental and Analytical Freezing Frac- 
tions Compared. n a calculated with length tl using 
equation (16) for Nu. 

±10% agreement band, although on average n e values 
tended to about 5% lower than the n a . The value of n e 
found from equation (1) is inversely proportional to the 
value of LWC through the accumulation parameter (eq. 
(2)). Because the accuracy for the LWC is estimated to 
be about ±12%, this agreement is as good as can be 
expected. 

While n a = 1 suggests fully rime ice, n e is susceptible to 
the uncertainties in both the LWC calibration and the 
tracing technique and can therefore be greater or less 
than unity. n e values for experimental rime shapes in 
this study were within ±10% of 1. 

Unlike the results of reference 5, the experimentally- 
determined freezing fractions for this study showed no 
systematic deviation from the analytical values at low 
freezing fractions. This consistently good agreement 
between n e and n a over a range of conditions including 
both fully rime and fairly warm glaze argues against the 
contention that the Messinger freezing fraction formula- 
tion neglects any important surface-water effects. 

Ice Shapes for Freezing Fractions Near 1 

For the conditions of 10-12-02, run 1 1, a static tempera- 
ture of -25°C was the maximum that would produce an 
analytical freezing fraction of 1. Four additional tests 
were made with temperatures of -22 to -18°C with cor- 
responding n a of 0.91 to 0.78 to determine the effect of 
freezing fraction on ice shape for conditions not quite 
fully rime. The results are shown in figure 4. The 
shape for fully rime ice is shaded in each portion of the 
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figure and the near-rime shape is given as a solid line 
for comparison. The rime shape conformed approxi- 
mately to the leading-edge contour of the model, as 
seen previously for somewhat different rime conditions 
in figure 2 (c). 

When the analytical freezing fraction was reduced from 
1 to 0.91 (fig. 4 (a)), a narrow strip of ice at the leading 
edge changed from the characteristic white appearance 
of rime to a slightly transparent form. At the same time 
the smooth, convex shape of typical rime within this 
narrow strip was replaced with a slight valley that can 
be seen in the figure. The remainder of the ice aft of 


this narrow strip at the leading edge maintained the 
rime appearance. Further reductions in freezing frac- 
tion to 0.87 (fig. 4 (b)), 0.83 (fig. 4 (c)) and 0.78 (fig. 4 
(d)) caused the leading-edge depression to deepen 
slightly without changing the rime appearance aft of the 
leading edge. 

These results show that the use of the Messinger analy- 
sis leads to a calculated freezing fraction of 1 when ice 
shapes are indeed fully rime and less than 1 when glaze 
features begin to show at the leading edge. Thus, the 
results of figures 3 and 4 together show that freezing 
fractions calculated according to equations (3) - (17) 



x/c 

n a =1.00, 10-12-02 Run 11 
n a = 0.91. 10-12-02 Run 8 

(a) Freezing Fraction of 0.91 Compared with Fully Rime 



x/c 

n a =1.00. 10-12-02 Run 11 
n a = 0.82. 10-11-02 Run 9 

(c) Freezing Fraction of 0.82 Compared with Fully Rime 



x/c 

n a =1.00, 10-12-02 Run 11 
n a =0.86, 10-12-02 Run 7 

(b) Freezing Fractions of 0.86 Compared with Fully 
Rime 



x/c 

w o =1.00. 10-12-02 Run 11 
n a =0.78. 10-1 1-02 Run 8 

(d) Freezing Fractions of 0.78 Compared with Fully 
Rime 


Figure 4. Effect of Freezing Fraction on Ice Shape at and near n = 1. See Table I for Test Conditions. 
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satisfy both tests (i) and (ii) of the Introduction. 

Effect of Modifications to the Freezing Fraction Calcu- 
lations 

Next the effect of modifying the freezing fraction calcu- 
lations will be shown. First, the comparison of experi- 
mental and analytical freezing fractions will be repeated 
with model chord substituted for d in equations (1), (2), 
(15) and (17). Note that the same value of n e is ob- 
tained from equation (1) whether c or d is used because 
both non-dimensional ice thickness and A c are inversely 
proportional to model size. Figure 5 (a) gives the re- 
sults of this exercise. The n a based on c were lower 
than those using d. Consequently, the mean n e was just 
over 10% higher than n a . The differences between n e 
and n a might still be explained by the L WC calibration 
uncertainty, so the general trend of the n e vs 1 n a com- 
parison does not indicate a problem. However for all 
three of the rime cases in Table I, replacing d with c 
resulted in n a decreasing to a value significantly below 
unity even though the experimental shapes obtained 
were fully rime. For example, n a for the conditions of 
figure 2(c), went from 1 to 0.82 and the rime reference 
conditions for fig. 4 now produced n a = 0.89. There- 
fore, the value of n a used in figure 5 cannot be correct if 
the LWC calibration is reasonably accurate. 

Figure 5 (b) shows another example of the effect of 
modifying the freezing fraction calculation, this time 
using d but calculating heat transfer from 

Nu = 0.0239Re°' 805 (21) 



Freezing Fraction from Analysis, n a 
(a) n u calculated with length c. 

Figure 5. Effect of Simple Modifications to Freezing 
Fraction Calculations 


instead of equation (16). Equation (21) was quoted by 
KreitlV from results published by Hilpert 13 for the aver- 
age Nu around a single cylinder in cross flow at Re in 
the range 40,000 to 400,000. For the present study Re 
varied from 55,000 to 120,000, so it might not seem 
unreasonable to use equation (21) to determine n a . 

Again, the values of n e were in reasonable agreement 
with this new value of n a as seen in figure 5 (b). How- 
ever, once again when test conditions produced rime 
shapes, this analysis produced freezing fractions below 
1. For example, the conditions of figure 2 (c) gave n a = 
0.89 and those of the rime reference shape for figure 4 
resulted in n a = 0.84. Thus, both the analyses used for 
figures 5a) and (b) failed test (ii) of the Introduction. 

The low value of n a at rime could be an indication of 
the inaccuracy of the tunnel L WC calibration. Suppose, 
for example, that the true tunnel LWC for the tests of 
Table I was 20% lower than the values reported there 
and used in the calculation of n a . Repeating the analy- 
sis of figure 5 (b) with all LWC’s reduced to 80% of the 
values in Table I, the values of n a increased about 15% 
such that rime conditions were now represented by an 
analytical freezing fraction of 1 or nearly 1. However, 
at the same time the n e were increased by about 25%, 
resulting in a much poorer comparison between the two 
freezing fractions, as can be seen in figure 5 (c). Fur- 
thermore, n e values for the conditions that produced 
rime shapes became roughly 15 - 20% greater than 1. 
This inconsistency between n a and n e indicates that the 
problems seen in figure 5 (b) are more fundamental 

s: 



Freezing Fraction from Analysis, n a 
(b) n a calculated using equation (21) with length d. 
Figure 5. (con’t). 
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Freezing Fraction from Analysis, n a 

(c) n a calculated with Nu used for figure 5 (b), but as- 
suming actual L WC 20% lower than calibration. 

Figure 5. (concluded) 

than a problem with the LWC calibration. 

Concluding Remarks 

Figure 3 illustrated that the two tests of the Introduction 
are satisfied by the Messinger analysis as represented 
by equations (3) - (17). (i) There was generally good 
agreement between the ice-thickness freezing fractions 
and the analytical values and (ii) the Messinger analysis 
led to freezing fractions of 1 for conditions that pro- 
duced rime shapes. Within the ability to measure ex- 
perimental ice shapes, there was no indication that the 
Messinger analysis overlooks any significant physical 
processes affecting the freezing fraction. 

The results in figure 5 show how leading-edge thick- 
ness measurements and the n e derived from them can be 
used to check both the soundness of the freezing frac- 
tion analysis and infer the reliability of the tunnel LWC 
calibration. 

These observations may not be valid for very short ic- 
ing sprays. For such situations, convection and conduc- 
tion heat transfer to and along the model surface may 
need to be included in the analytical expression for 
freezing fraction. Some analysis is needed to estimate 
how far into the icing spray such terms might be sig- 
nificant. 


For scaling studies, scale and reference values of freez- 
ing fraction are matched to establish scale temperature 
or LWC. As long as the same analytical expressions are 
used for both reference and scale, some variation in the 
freezing fraction analysis is possible with little effect on 
the calculated scale conditions. For applications such 
as ice-accretion prediction codes, it is more important 
to insure the correct analytical expressions. 

At freezing fractions of 0.3 and less, variations in re- 
corded leading-edge thickness from one span-wise loca- 
tion to another of as much as 30% were observed. At 
higher freezing fractions the variations were typically 
less than 8%. Much of the difference noted in ice 
thickness was an indication of the error inherent in 
manual tracing of ice. To some extent, the random na- 
ture of ice accretion is also revealed by these tracings. 
Currently used tracing methods need to be improved for 
better analysis, although these methods are adequate for 
routine ice-shape comparisons and scaling studies. 

For this work, the Messinger analysis for freezing frac- 
tion as described by Ruff was applied using twice the 
leading-edge airfoil radius as the relevant length for the 
various parameters. Results using chord produced ana- 
lytical freezing fractions as low as 0.82 for conditions 
that gave fully rime ice shapes. This inconsistency 
suggested that the use of chord for the relevant length is 
not correct. 

A convective heat-transfer film coefficient represented 
by Nu Re 0> was used in the analysis. Analytical 
freezing fractions determined using an alternate form, 
Nu Re o s , also resulted in freezing fractions of less 
than 1 at rime conditions. 

For the Appendix-C conditions at which this study was 
performed, the 1RT LWC calibration was consistent 
with experimental ice shapes and freezing fractions. 

References 

1 Messinger, B.L., “Equilibrium Temperature of an Un- 
heated Icing Surface as a Function of Airspeed,” J. 
Aeron. Sci. vol. 20 no. 1, January 1953, pp 29-42. 

2 Ruff, G.A., “Analysis and Verification of the Icing 
Scaling Equations,” AEDC-TR-85-30, vol 1 (Rev), 
March 1986. 

3 Anderson, David N., “Acceptable Tolerances for 
Matching Icing Similarity Parameters in Scaling Appli- 
cations,” A1AA-200 1-0832, January 2001. 

4 Langmuir, Irving and Blodgett, Katharine B.: “A 
Mathematical Investigation of Water Droplet Trajecto- 
ries,” Army Air Forces Technical Report No. 5418, 
February 1946. 


NAS A/CR— 2005-2 13852 


10 




5 Bilanin, A.J., “Proposed Modifications to the Ice Ac- 
cretion/Icing Scaling Theory,” AIAA Paper 
AIAA-88-0203, January 1988. 

6 Tribus, Myron, Young, G.B.W. and Boelter, L.M.K., 
“Analysis of Heat Transfer Over a Small Cylinder in 
Icing Conditions on Mount Washington,” Trans. ASME 
vol. 70, November 1948, pp 971-976. 

7 Charpin, Francois and Fasso, Guy, “Essais de givrage 
dans la grande soufflerie de Modane sur maquettes a 
echelle grandeur et echelle reduite,” L ’Aeronautique et 
I’Astronautique, no. 38, 1972, pp 23 - 31. English 
translation published as “Icing Testing in the Large 
Modane Wind-Tunnel on Full-Scale and Reduced Scale 
Models,” NASA TM-75373, March 1979. 

8 Kreith, Frank, Principles of Heat Transfer, Interna- 
tional Textbook Co., Scranton, 1958. 


9 Poinsatte, Philip E., “Heat Transfer Measurements 
from a NACA 0012 Airfoil in Flight and in the NASA 
Lewis Icing Research Tunnel,” NASA CR 4278, March 
1990. 

10 Ide, Robert F. and Oldenburg, John R., “Icing Cloud 
Calibration of the NASA Glenn Icing Research Tun- 
nel,” A1AA-200 1-0234, January 2001. 

11 Anderson, David N., “A Preliminary Study of Ice- 
Accretion Scaling for SLD Conditions,” AIAA-2002- 
0521, January 2002. 

12 Anderson, David N., and Feo, Alejandro “Ice- 
Accretion Scaling Using Water-Film Thickness Pa- 
rameters,” A1AA-2002-0522, January 2002. 

13 Hilpert, R., “Warmeabgabe von geheizten Drahten 
und Rohren,” Forsch. Gebiete Ingenieurw., vol. 4, 
1933, p 215. 


NASA/CR- 


-2005-213852 


11 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1 . AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

August 2005 

3. REPORT TYPE AND DATES COVERED 

Final Contractor Report 


4. TITLE AND SUBTITLE 


Evaluation and Validation of the Messinger Freezing Fraction 


6. AUTHOR(S) 

David N. Anderson and Jen-Ching Tsao 


5. FUNDING NUMBERS 


WBS-22-728^11-17 

NCC3-884 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Ohio Aerospace Institute 
22800 Cedar Point Road 
Brook Park, Ohio 44142 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 15222 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA CR— 2005-213852 
A1AA-2003-1218 


11. SUPPLEMENTARY NOTES 

Prepared for the 41st Aerospace Sciences Meeting and Exhibit sponsored by the American Institute of Aeronautics and 
Astronautics, Reno, Nevada, January 6-9, 2003. Project Manager, Thomas H. Bond, Instrumentation and Controls 
Division, NASA Glenn Research Center, organization code RI, 2 1 6-433-3900. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category: 03 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390. 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

One of the most important non-dimensional parameters used in ice-accretion modeling and scaling studies is the freezing 
fraction defined by the heat-balance analysis of Messinger. For fifty years this parameter has been used to indicate how 
rapidly freezing takes place when super-cooled water strikes a solid body. The value ranges from 0 (no freezing) to 1 
(water freezes immediately on impact), and the magnitude has been shown to play a major role in determining the physical 
appearance of the accreted ice. Because of its importance to ice shape, this parameter and the physics underlying the 
expressions used to calculate it have been questioned from time to time. Until now, there has been no strong evidence either 
validating or casting doubt on the current expressions. This paper presents experimental measurements of the leading-edge 
thickness of a number of ice shapes for a variety of test conditions with nominal freezing fractions from 0.3 to 1.0. From 
these thickness measurements, experimental freezing fractions were calculated and compared with values found from the 
Messinger analysis as applied by Ruff. Within the experimental uncertainty of measuring the leading-edge thickness, 
agreement of the experimental and analytical freezing fraction was very good. It is also shown that values of analytical 
freezing fraction were entirely consistent with observed ice shapes at and near rime conditions: At an analytical freezing 
fraction of unity, experimental ice shapes displayed the classic rime shape, while for conditions producing analytical 
freezing fractions slightly lower than unity, glaze features started to appear. 


14. SUBJECT TERMS 

Icing physics; Scaling; SLD; LWC 


15. NUMBER OF PAGES 

17 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 




